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Abstract 

The caloric curve (excitation energy per particle as a function of temperature) 
for finite nuclei is calculated within the non-linear Walecka model for different 
proton fractions and different parameterizations. The results obtained are compared 
with published experimental data and other theoretical results. Other properties 
related with the droplet formation as the surface energy, neutron skin thickness 
and binding energy per nucleon are also displayed as a function of temperature and 
their behaviours are discussed. It is shown that the caloric curve is sensitive to the 
proton fraction and to the Coulomb interaction. At T = the droplet properties 
do not depend on the parameterization used. At finite temperature however, the 
sensitivity appears. 



PACS number(s): 21.10.-k, 21.30.-x,21.65.+f,25.70.-z 
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1 Introduction 



One of the most important problems in contemporary nuclear physics and astrophysics is 
the determination of the properties of nuclear matter as functions of density, temperature 
and the neutron-proton composition. This information is important in understanding 
the explosion mechanism of supernova, the cooling rate of neutron stars or the processes 
involved in the formation of the trans-iron elements |1|, |2|, |3|, |], [5|. |(| . 

Recent advances in experiments using heavy-ions at high energies and radioactive 
beams with large neutron or proton excess have made it possible to create not only nuclei 
at the limits of stability, but also a transient state of nuclear matter with appreciable 
isospin asymmetry, thermal excitation, and compression [0, || [|. 

The production of several intermediate mass fragments in a short time scale dur- 
ing heavy ion collisions is known as nuclear multifragmentation. In multifragmentation 
experiments, an equilibrated system is always formed and the behaviour of the energy 
fluctuations suggests that the system undergoes a liquid-gas phase transition. One of 
the evidences of this transition is the fact that the heat-capacity exhibits a peak at a 
certain temperature. Although different pictures are normally used in order to describe 



this phase transition phenomenon [JTOj, [11], |12j] , some common features are present in all of 
them. The caloric equation of state, which is given by the excitation energy per nucleon in 
terms of the thermodynamic temperature is an important quantity to be investigated in 
the search for a phase transition. Notice that the temperature at which phase transitions 
take place experimentally is not the same as the thermodynamic one. 

In order to explain recent experimental results, one must study not only the ground 
and excited states of normal nuclei, but also nuclear states of high excitation and far 
from stability. In particular it is important to understand the role of the isospin degree 
of freedom in heavy-ion collisions at intermediate energies and the properties of nuclear 
matter in the region between symmetric nuclear and pure neutron matters. 

The radius, thickness of neutron skins, deformation, binding energy, density distribu- 
tions, and other properties of radioactive nuclei near the drip lines depend sensitively on 
the isospin-dependence of the nuclear equation of state (eos). The authors of references 
H| and |14| have recently stressed the possibility of extracting the EOS of asymmetric 



nuclear matter through the investigation of these properties as well as via reactions of 
neutron-rich nuclei at intermediate energies. 
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The properties of neutron rich nuclei have been studied using, among others, Skyrme- 
Hartree-Fock (shf) and the relativistic mean- field (rmf) theories [|T5[]. In particular it 



has been shown by Tanihata |13[ that in shf the saturation density essentially does 
not change, but in rmf it decreases rapidly as nuclear matter becomes more neutron- 
rich. Most EOS used in the description of neutron stars and supernova are obtained 
from nonrelativistic models. It is therefore of interest to study the EOS, namely the 
inhomogeneous phase at the liquid-gas phase transition, using relativistic models. The 
liquid-gas phase transition is probed not only by investigating the properties of neutron 
stars but also through experimental results obtained in collisions involving heavy-ions and 
in multifragmentation experiments, already mentioned. 

Nucleation process, which is known as the liquid droplet formation in a gas back- 
ground, is a very old problem with many applications in different areas of physics and 
chemistry. Newton (1687) and Laplace (1816) dedicated themselves to understand the 
physics underlying the nucleation theory earlier than the formulation of the zeroth law 
of thermodynamics in 1913 [[HJ. In this work we turn into the nucleation process inside 
nuclear matter. 

Within the framework of relativistic models, the liquid-gas phase transition in nuclear 
matter has been investigated at zero and finite temperatures for symmetric and asym- 
metric semi- infinite systems [17] , |18| , I9| , |20| , |2l|j . Droplet formation in the liquid-gas phase 
transition in cold [23] and hot asymmetric nuclear matter in the context of rela- 
tivistic mean field theory, namely using the non-linear Walecka model (NLWM) p6f has 



been studied recently using a Thomas Fermi approximation. As shown in Refs. 24 



the optimal nuclear size of a droplet in a neutron gas is determined by a delicate balance 
between nuclear Coulomb and surface energies. The surface energy favors nuclei with a 
large number of nucleons A, while the nuclear Coulomb self-energy favors small nuclei. 



In a previous work [^5[], the influence of the proton fraction in the caloric curve was 
calculated for nuclei obtained within the framework of the approach mentioned above. 
In studying multifragmentation, an input parameter called the freeze-out volume, which 



simulates a phase transition at constant volume, is normally used fljj| . The consequences 
for the caloric curve of imposing thermalization in a freeze-out volume have also been 
discussed. 

The present work extends our study, by investigating the influence of the temperature 
in some important quantities related to the droplet formation, such as the surface energy, 
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the neutron skin thickness and the binding energy per nucleon. The importance of the 
Coulomb interaction and the nucleus proton fraction is investigated. The effects of using 
different parameterizations are also analysed. In this way, a more complete understanding 
of the caloric curve is obtained. 

The paper is organized as follows. In section 2 the thermodynamical potential in the 
framework of the Thomas-Fermi approximation is calculated for the NLWM and the two- 
phase coexistence is discussed. In section 3 we present the numerical results. Finally, in 
the last section the conclusions are drawn. 



2 The Thomas- Fermi Approximation in the Extended 
Non Linear Walecka Model 



In what follows we describe the equation of state of asymmetric matter within the frame- 
work of the relativistic non-linear Walecka model [2£| , [27]] with the inclusion of p- mesons 
and the electromagnetic field. In this model the nucleons are coupled to scalar-isoscalar 
0, vector-isoscalar V^, vector-isovector meson fields and the electromagnetic field A 11 . 
The lagrangian density reads: 



£ = ip 



eN 1 



{M-g s 



1> 



1 

I 



1 

_< 

4' 



1 
2 



where Q 



dubf, 



\mlb„ ■ - -F^F^ 



9 P {bp x K) and 



d v A p . The 



model comprises the following parameters: three coupling constants g s , g v and g p of the 
mesons to the nucleons, the nucleon mass M, the masses of the mesons m s , m v , m p , the 
electromagnetic coupling constant e = \J~^§f and the self-interacting coupling constants k, 
A and £. In this work, unless otherwise stated, we use a set of constants usually identified 
as NL1 [p8|| , which is not unique, but it gives a good description of the ground state 
properties of many stable nuclei. The consequences of using different parameterizations 
is also investigated and, for this purpose, we use two other sets known as NL3 and 
TM1 [30|. The values of the constants are displayed in table 1 and the most important 
bulk properties obtained with these three parameterizations are displayed in table la. 
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The basic quantity in the Thomas-Fermi approximation is the phase-space distribu 
tion function for protons and neutrons: 

1 



/i±(r,p,t) 



, i = p,n, 



(2) 



1 + exp[(e =F Vi)/T] 

where Vi — ^ — Vjo are the effective chemical potentials with /ij being the chemical 
potentials for particles of type i, 



9p 



9 Pl 



V p o = 9vV + —b + eA , V„ = g v V - —b , 

e = \/ p 1 + M* 2 , M* = M — g s (p is the effective nucleon mass and T is the temperature. 
The classical entropy of a Fermi gas is given by 



s - -2 E / {f»* (t^w) + ^ ~ m + (/. + « 

and the thermodynamic potential is defined as 



(3) 



fl = E - TS - Vi N i- 



i=p,n 



where N p , N n are, respectively, the proton and the neutron number: 



N; 



d 3 rp*M), Pl = 2 
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From the above expressions we get for ([|), in the static approximation, 

O = J d 3 r Q [(V</>) 2 - (W ) 2 - (V&o) 2 - (VA ) 2 ] - ^ 
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The energy density reads 
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For later discussion it is convenient to separate £(r,T) into two parts: the nuclear con- 
tribution, £ N , and the electromagnetic one; the latter is 



^coul. 



^eA (/ p+ -/ p _)-i(VA ) 2 



The fields that minimize Q satisfy the equations 

9V ef 
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where p B = p P + p n , Ps = P P ~ Pn and 

i=p,n ' 



d 3 p M* 
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(12) 



These coupled differential equations are solved numerically and all relevant quanti- 
ties which depend on the fields are calculated. Looking for phase transitions in binary 
systems such as this one requires the study of three kinds of instabilities which can oc- 
cur: mechanical, diffusive and thermodynamical. The condition for mechanical stability 
requires 

dP 



>0, 



(13) 



PPbj y p ,t 

where P is the pressure and Y p = p p / Pb is the proton fraction. The condition for diffusive 
stability implies the inequalities 

'dp p \ ^ n ( dpr, 



> and 
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These conditions reflect the fact that in a stable system, energy is required to increase the 
proton concentration at constant pressure and temperature. Thermo dynamical stability 
is expressed by 

*=(£L*°" (i5) 

where c v is the specific heat and 

e* = e(T) -e(T = 0), (16) 

is the excitation energy per particle. The total energy per particle at temperature T is 
given by El 



e(T) = J ^ld 3 r = e N (T) + e coul .(T) , (17) 

where A is the total number of particles in the volume under consideration, defined at the 
end of section 3, and En(T), e cou l(T) denote, respectively, the nuclear and electromagnetic 
contributions. 

The two-phase liquid-gas coexistence is governed by the Gibbs condition 

fii(p p , p n , M*) = tii(p' p ,p' n ,M*'), i=p,n (18) 
P(p p ,p n ,M*) = P(p' p , p' n , M*') , (19) 

where the primed and unprimed quantities correspond to the two different phases. 

In the mean field approximation for infinite nuclear matter, the meson fields are re- 
placed by their expectation values. From the equations of motion, they can easily be 
obtained [^] and the thermodynamic quantities of interest are given in terms of these 
meson fields. We have made use of the geometrical construction |J2], in order to obtain 
the chemical potentials in the two coexisting phases for each pressure of interest. In order 
to obtain the binodal section which contains points under the same pressure for different 
proton fractions, we have solved eqs. (|18| ) and (|19| ) simultaneously with the following 
ones: 

k X 

m 2 s (j) + -(j) 2 + -(f)l = g s p s (iy p ,u n ,M*) (20) 



and 



m 2 M + + ^0o' 3 = 9s P s{v' p , K, M*'). (21) 
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3 Numerical Results for Finite Systems 



The solution for the infinite system gives us the initial and boundary conditions for the 
program which integrates the set of coupled non-linear differential equations ([D to (^) 
in the Thomas-Fermi approximation. In this work the numerical calculation was carried 
out with the iteration procedure described in Refs. W\ [23], |2l . 



Some quantities of interest are outlined below. The surface energy per unit area of 
the droplets in the small thickness approximation is 



a 



dr 



d(f) 
dr 



dVo 
dr 



dr 



The proton and neutron radii in the spherical geometry, Ri 

r# 1 



pi{r)r dr 



(22) 

i — p,n), are defined as 
(23) 



where p i; i and p i)9 refer to the liquid and gas density respectively and R' is the value of r 
for which the fields and density reach their asymptotic gas values. 
The neutron skin thickness is defined as [E0| 



Q = R ri 



R p . 



(24) 



Other quantities of interest are the number of protons (Z) and neutrons (N) in the droplet, 
given by 



and 



47T 



N = Air 



R' 



p p (r)r 2 dr - ^-R' 3 p p<g 



R> 



p n (r)r 2 dr - ^i? /3 pn, 9 - 



(25) 



(26) 



e(T) - M, 



(27) 



The binding energy per nucleon is 

B 
A 

where A = Z + N. 

In the sequel we mainly study two nuclei: g2°Sm 88 and g2 6 Sm 10 4. We have chosen these 
two nuclei because they lie in the mass range of interest for intermediate energy heavy- 
ion collisions [0, On the other hand, two isotopes with quite different number of 
neutrons are taken in order to study the effect of proton-neutron asymmetry. 
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In figure [11 the proton density profiles are plotted with and without the electromagnetic 
field for the g2°Sm 88 at T = 0. The central density is smaller when the Coulomb interaction 
is considered because more protons are pushed off to the surface. The same behaviour 
is obtained for any other fixed temperature, whenever both systems, with and without 
the inclusion of the electromagnetic field are found. Notice, however, that the critical 
temperature, from which the droplet ceases to exist, is smaller when the electromagnetic 
field is included, as compared with systems where the Coulomb interaction is disregarded. 

In figure |2|, the g| Smgs density profiles for the protons are plotted with the inclusion 
of the electromagnetic field, for temperatures varying from T = to T = 6.5 MeV. The 
central density decreases with the increase of the temperature and consequently, the hotter 
the droplet, the larger it is. The same behaviour is observed when the Coulomb interaction 
is not included and for any other proton fraction. These results are in accordance with 



those obtained in pi |, with a non-relativistic model. 

In table 2, we show the surface energy, the proton radius, the neutron skin thickness, 
the central density, the binding energy per nucleon and the excitation energy per particle 
calculated without the inclusion of the Coulomb interaction for the g2°Sm 88 , which has 
a proton fraction equal to 0.41. In table 3, the same quantities are shown when calcu- 
lated with the electromagnetic interaction included. In tables 4 and 5, these quantities 
are obtained for the Q2 6 Sniio4, which has a proton fraction of 0.37, without and with the 
Coulomb interaction respectively. Examining these tables, several common trends are 
noticed, despite the difference in the proton fraction. In general, the surface energy, the 
central density and the binding energy decrease with the increase of the temperature. 
Inversely, the proton radius and the excitation energy per particle increase as the temper- 
ature becomes higher. Comparing tables 2 with 3 and 4 with 5, one observes that the non 
inclusion of the Coulomb interaction gives rise to smaller and denser nuclei with a bigger 
neutron skin: the Coulomb force pushes the protons to the surface as already mentioned. 
It is also seen that the binding energy is always larger for the systems considered with- 
out the Coulomb interaction, in the sense that a more negative binding energy reflects a 
more bound state. This fact is easily explained since the Coulomb interaction tends to 
move protons apart and hence, diminishes the binding energy. It is however, interesting 
to notice that the excitation energy increases faster with temperature when the Coulomb 
energy is taken into account. Being a long range force one could think that there would 
be essentially no influence on the excitation energy. In fact, for the full calculation, one 
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has £*{T) = e(T) — e(0) ~ £n(T) — £jv(0), which is explained by the long range of the 
Coulomb force. If the calculation is performed with no Coulomb force (NC), the exci- 
tation energy e* NC (T) = e^ c (T) — e^ c (0) contains only the nuclear contribution. In 
figure ^ the effect of the electromagnetic force is clear: if the Coulomb field is taken into 
account the excitation energy is higher, independently of the proton-neutron asymmetry, 
i.e. e* NC (T) < e*(T). Indeed, we should not forget that the Coulomb force affects the 
proton distribution function and, therefore, all properties of the nucleus. As already men- 
tioned before, we remind that in tables 3 and 5, the last temperature (6.5 MeV) is the 
one at which the droplets cease to exist within our framework, since we are not able to 
obtain convergence for a droplet of the size considered at higher temperatures; the same 
is not true in tables 2 and 4, where droplets can still be found for temperatures larger 
than 9 MeV. 

In tables 6 and 7, we compare for the nuclei Q2°Sm 88 and Q2 6 Smi 4, respectively, the 
results obtained for the quantities already mentioned above, when different parameteri- 
zations are used. We have considered three values of the temperature and included the 
Coulomb interaction. As previously noticed for the parameterization NL1, also for the 
other two parameterizations, the surface energy, the central density and the binding en- 
ergy per particle decrease with the increase of the temperature, while the proton radius 
and the excitation energy per particle increase, under these conditions. 

At T = 0, the results given by all the parameterizations considered are quite similar. 
This reflects the fact that all the parameterizations have been fitted to the ground-state 
properties of stable (NL1) or stable and unstable (NL3 and TM1) nuclei. However, for 
finite temperatures some differences must be referred. Although at T = the surface 
energy is essentially the same for the 3 parameterizations, it is seen that for NL1 it 
reduces much faster with T, while in the TM1 it suffers the smallest reduction. This is 
true for both proton fractions, Y p = 0.41 and 0.37. In the last case, however, a reduces 
faster with T for all parameterizations. 

The neutron skin thickness, G, is not sensitive to the temperature, but it is interesting 
to observe that, for a given T, it assumes smaller values for NL3 and TM1 than for NL1. 



This behaviour was already discussed in [g7|] for T = and is related with the smaller 
values predicted to the symmetry energy for nuclear matter in NL3 (37.4 MeV) and in 
TM1 (36.9 MeV) as compared to NL1 (43.5 MeV). Also related with these values are the 
slightly higher binding energies and the slowlier increasing of the excitation energies with 
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temperature for the NL3 and TM1 parameterizations. 

In tables 8 and 9 several properties of isotopes with Z = 62 are listed for T = 
and T = 5 MeV, respectively. One conclusion can immediately be taken: the proton- 
neutron asymmetry makes the nuclei softer at higher temperatures and properties like 
the surface energy, the central density or the proton distribution radius vary much more 
with the number of neutrons if temperature is different from zero. It is important that 
a parameterization of the surface energy with temperature takes into account the isospin 
dependence. The decrease of the central density with the decrease of the proton fraction 
is typical of relativistic mean-field models. This feature is different from the prediction 



of the Skyrme-Hartree-Fock calculations [|14]], where the central density does not change 
with the proton fraction. We confirm that temperature has no effect on the neutron skin. 
In figure we show the proton density profiles for Z = 62 and different proton fractions, 
at T = 5 MeV. We observe that the proton central density increases with the proton 
fraction. 

In figure [5] we show the caloric curves obtained with the Coulomb interaction, for 
g2°Sm 88 (thick full line) and g2 6 Sm 104 (thick dashed line), the excitation energies at 
T = 5 MeV, for different proton fractions (big triangles, from the left to the right 
Y p = 0.33,0.35,0.37,0.39,0.41 and 0.44), the experimental data from |7j (diamonds) and 
[ |35|j (stars), and the Fermi-gas law e* = 1/kT 2 , with k = 13.0 (thin dashed line). We 
conclude that the excitation energy for 166 Sm (thick dashed curve), proton fraction 0.37, 
increases slowlier with temperature than for 150 Sm (thick full curve), proton fraction 
0.41. These two curves are consistent with data of [|35] and a level density parameter 
A/k, k = 13.0 in the Fermi gas model relation. The calculation at T = 5 MeV, for the 
proton fractions represented in the figure, shows that the caloric curve is sensitive to the 
proton-neutron ratio in the compound nucleus. 



4 Conclusions 

In summary, we have studied the surface properties and the excitation energies of arising 
droplets in a vapor system for temperatures up to 6.5 MeV. The droplets are described 
in terms of a Walecka-type model within the Thomas-Fermi approximation. It was 
shown that, for a fixed temperature, nuclei have a softer behaviour with the decrease 
of its proton fraction, namely the surface energy and the central density are smaller 
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and the proton radius becomes larger. We have also concluded that the neutron skin 
is independent of the temperature. Another important conclusion refers to the effect of 
the Coulomb interaction: it reduces the central density of nuclei, as well as the neutron 
skin, since the protons are pushed out to the surface, and increases the excitation energy. 
We have also compared the behaviour with temperature of the results obtained within 
three different parameterizations of the non-linear Walecka model, as previously referred, 
NL1 [p8[| , NL3 ]2l| , TM1 [30|. It was shown that, at T = 0, the three parameterizations 
give similar results but, at finite temperature, some of the properties investigated exhibit 
different behaviours: in the NL1 parameterization, the nuclei reacts in a softer way to 
temperature, while in the TM1 parameterization the properties of nuclei change slowlier 
with temperature. The excitation energies of droplets either corresponding to 150 Sm or 
166 Sm, for temperatures between 3 and 6.5 MeV, are consistent with the caloric curve in 
the Fermi gas approximation with a level density parameter A/13. This result agrees with 



experimental data obtained in heavy-ion collisions at intermediate energies |35| . We show 
that the caloric curve is sensitive to the proton fraction and therefore to the symmetry 
term of the model used. It is also shown that close to the critical temperature the three 
parameterizations tested give different results. Experimentally the dependence on the 
proton fraction could be studied by comparing data obtained from sources with different 
proton fractions. 
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Table 1. Sets of parameters used in this work. All masses are given in MeV. 



Force 


[Ref.] 


M 


m s 


m v 


m p 


9s 


9v 


9 P 


k/M 


A 




NL1 


[28] 


938. 


492.25 


795.36 


763.0 


10.138 


13.285 


9.952 


5.122 


-217.613 


0.0 


NL3 


[29] 


939. 


508.194 


782.501 


763.0 


10.217 


12.868 


8.948 


4.377 


-173.31 


0.0 


TM1 


[30] 


938. 


511.198 


783.0 


770.0 


10.0289 


12.6139 


9.2644 


3.04 


3.7098 


0.0169 



Table la Nuclear matter properties in the context of several parameterizations. 





NL1 [28| 


NL3 [29] 


TM1 [30] 


B/A (MeV) 


16.5 


16.3 


16.3 


Po (fm~ 3 ) 


0.153 


0.148 


0.145 


K (MeV) 


211 


272 


281 


a sym , (MeV) 


43.7 


37.4 


36.9 


M*/M 


0.57 


0.60 


0.63 



Table 2. Output results given by the solution of the coupled differential equations for 
different temperatures without the inclusion of the electromagnetic field for g^Snigs (Y p = 
0.41) with the NL1 parameterization. 



T 


a 


R p 


e 


P(0) 


B/A 


e*(T) 


(MeV) 


(MeV fm" 2 ) 


(fm) 


(fm) 


(fm" 3 ) 


(MeV/ A) 


(MeV/A) 


0. 


1.15 


5.88 


0.33 


0.161 


-11.82 


0.00 


3. 


1.14 


5.77 


0.39 


0.166 


-10.95 


0.87 


4. 


1.06 


5.84 


0.37 


0.163 


-10.64 


1.18 


5. 


0.95 


5.92 


0.36 


0.156 


-10.14 


1.68 


6. 


0.83 


6.04 


0.35 


0.147 


-9.49 


2.33 


6.5 


0.74 


6.13 


0.35 


0.140 


-8.99 


2.83 


7. 


0.68 


6.22 


0.35 


0.137 


-8.59 


3.23 


8. 


0.54 


6.41 


0.34 


0.126 


-7.56 


4.26 


9. 


0.40 


6.63 


0.36 


0.112 


-6.33 


5.49 
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Table 3. Output results given by the solution of the coupled differential equations for 
different temperatures with the inclusion of the electromagnetic field for g2°Sm 88 (Y p = 
0.41) with the NL1 parameterization. 



T 


a 


R p 





P(0) 


B/A 


e*(T) 


(MeV) 


(MeV fm" 2 ) 


(fm) 


(fm) 


(fm~ 3 


(MeV/A) 


(MeV/A) 


0. 


1.00 


6.22 


0.16 


0.142 


-8.39 


0.00 


3. 


0.98 


6.23 


0.18 


0.140 


-7.56 


0.83 


4. 


0.91 


6.30 


0.17 


0.137 


-7.06 


1.33 


5. 


0.81 


6.40 


0.16 


0.131 


-6.37 


2.02 


6. 


0.65 


6.57 


0.15 


0.121 


-5.38 


3.01 


6.5 


0.57 


6.68 


0.14 


0.115 


-4.78 


3.61 



Table 4. Output results given by the solution of the coupled differential equations for 
different temperatures without the inclusion of the electromagnetic field for g2 6 Sm 104 (Y p = 
0.37) with the NL1 parameterization. 



T 


a 


R p 


O 


P(0) 


B/A 


e*{T) 


(MeV) 


(MeV fm' 2 ) 


(fm) 


(fm) 


(fm" 3 ) 


(MeV/A) 


(MeV/A) 


0. 


1.03 


6.02 


0.51 


0.157 


-11.03 


0.00 


2. 


0.99 


6.03 


0.53 


0.154 


-10.71 


0.32 


3. 


0.92 


6.08 


0.54 


0.149 


-10.43 


0.60 


4. 


0.84 


6.16 


0.54 


0.144 


-10.12 


0.91 


5. 


0.74 


6.25 


0.52 


0.138 


-9.70 


1.33 


6. 


0.61 


6.39 


0.53 


0.129 


-9.12 


1.91 


7. 


0.49 


6.57 


0.53 


0.120 


-8.44 


2.59 


8. 


0.35 


6.87 


0.55 


0.105 


-7.57 


3.46 


9. 


0.27 


7.08 


0.52 


0.098 


-6.46 


4.57 
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Table 5. Output results given by the solution of the coupled differential equations for 
different temperatures with the inclusion of the electromagnetic field for g2 6 Smi 4 (Y p = 
0.37) with the NL1 parameterization. 



T 


a 


R p 





P(0) 


B/A 


e*(T) 


(MeV) 


(MeV fm" 2 ) 


(fm) 


(fm) 


(fm" 3 ) 


(MeV/ A) 


(MeV/A) 


0. 


0.92 


6.40 


0.32 


0.138 


-7.93 


0.00 


2. 


0.91 


6.42 


0.31 


0.137 


-7.65 


0.28 


3. 


0.86 


6.46 


0.32 


0.134 


-7.26 


0.67 


4. 


0.77 


6.54 


0.33 


0.129 


-6.80 


1.13 


5. 


0.66 


6.67 


0.32 


0.122 


-6.14 


1.79 


6. 


0.53 


6.87 


0.31 


0.113 


-5.20 


2.73 


6.5 


0.47 


7.00 


0.29 


0.108 


-4.81 


3.12 
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Table 6. Output results given by the solution of the coupled differential equations with the inclusion of the electromagnetic 
field, for Q2°Sm 88 (1^=0.41). For different temperatures, the results obtained with several parameterizations are shown. 







T=0 MeV 






T=5 MeV 






T=6 MeV 






NL1 


NL3 


TM1 


NL1 


NL3 


TM1 


NL1 


NL3 


TM1 


a (MeV fin" 2 ) 


1.000 


1.007 


0.987 


0.810 


0.857 


0.874 


0.650 


0.725 


0.762 


R P (fm) 


6.22 


6.26 


6.30 


6.40 


6.37 


6.39 


6.57 


6.49 


6.50 


(fm) 


0.16 


0.11 


0.09 


0.16 


0.12 


0.10 


0.15 


0.11 


0.10 


p(0) (fm" 3 ) 


0.142 


0.142 


0.139 


0.131 


0.134 


0.133 


0.121 


0.127 


0.127 


B/A (MeV/ A) 


-8.39 


-8.43 


-8.52 


-6.37 


-6.49 


-6.66 


-5.38 


-5.64 


-5.85 


e*(T) (MeV/ A) 


0.00 


0.00 


0.00 


2.02 


1.93 


1.87 


3.01 


2.79 


2.67 



Table 7. Output results given by the solution of the coupled differential equations with the inclusion of the electromagnetic 
field, for g2 6 Sniio4 (Y p =0.37). For different temperatures, the results obtained with several parameterizations are shown. 







T=0 MeV 






T=5 MeV 






T=6 MeV 






NL1 


NL3 


TM1 


NL1 


NL3 


TM1 


NL1 


NL3 


TM1 


a (MeV fin" 2 ) 


0.920 


0.933 


0.917 


0.660 


0.723 


0.742 


0.530 


0.615 


0.640 


R P (fm) 


6.40 


6.43 


6.49 


6.67 


6.62 


6.66 


6.87 


6.75 


6.76 


(fm) 


0.32 


0.25 


0.22 


0.32 


0.26 


0.23 


0.31 


0.25 


0.22 


p(0) (fm" 3 ) 


0.138 


0.139 


0.137 


0.122 


0.128 


0.127 


0.113 


0.120 


0.122 


B/A (MeV/ A) 


-7.93 


-8.04 


-8.11 


-6.14 


-6.30 


-6.43 


-5.20 


-5.50 


-5.67 



e*(T) (MeV/A) 0.00 0.00 0.00 1.79 1.74 1.68 2.73 2.54 2.45 



Table 8. Output results given by the solution of the coupled differential equations for 
T = with the inclusion of the electromagnetic field for Z = 62 and different number of 
neutrons with the NL1 parameterization. 



N 


Y 


a 


R p 


6 


P(0) 


B/A 






(MeV far 2 ) 


(fm) 


(fm) 


(fm" 3 ) 


(MeV/A) 


78 


0.445 


1.05 


6.11 


0.05 


0.144 


-8.59 


88 


0.413 


1.00 


6.22 


0.16 


0.142 


-8.39 


94 


0.397 


0.97 


6.29 


0.22 


0.141 


-8.25 


104 


0.374 


0.92 


6.40 


0.32 


0.138 


-7.93 


112 


0.358 


0.87 


6.50 


0.39 


0.136 


-7.76 


124 


0.334 


0.80 


6.62 


0.50 


0.133 


-7.37 



Table 9. Output results given by the solution of the coupled differential equations for 
T = 5 with the inclusion of the electromagnetic field for Z = 62 and different number of 
neutrons with the NL1 parameterization. 



N 


Y 


a 


R p 





P(0) 


B/A 


e*(T) 






(MeV fm" 2 ) 


(fm) 


(fm) 


(fm" 3 ) 


(MeV/A) 


(MeV/ A) 


78 


0.445 


0.87 


6.26 


0.04 


0.134 


-6.32 


2.27 


88 


0.413 


0.81 


6.40 


0.16 


0.131 


-6.37 


2.02 


94 


0.396 


0.75 


6.46 


0.23 


0.128 


-6.32 


1.93 


104 


0.374 


0.66 


6.67 


0.32 


0.122 


-6.14 


1.79 


112 


0.355 


0.55 


6.89 


0.39 


0.115 


-6.13 


1.63 


124 


0.334 


0.41 


7.27 


0.48 


0.102 


-6.01 


1.37 
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Figure 1: The g2°Sm 88 (Y p = 0.41) density profiles for protons p p (r) in fm 3 are plotted 
with (solid) and without (dashed) the inclusion of the electromagnetic field for T = 
MeV. 
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Figure 2: The g2°Sm 88 (Y p = 0.41) density profiles for the protons p p (r) in fm~ 3 are plotted 
with the inclusion of the electromagnetic field for temperatures varying from T = (at 
the top, on the left) to T = 6.5 MeV (at the bottom, on the left). 
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Figure 3: The caloric curves are shown for g2°Sm 88 (Y p = 0.41) with (thick solid line) and 
without (thin solid line) the Coulomb interaction and for g2 6 Sm 10 4 (Y p = 0.37) with (thick 
dashed line) and without (thin dashed line) the Coulomb interaction. 
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Figure 4: Density profiles for the protons p p {r) in fm 3 are plotted for T = 5 MeV 
and Y p = 0.44 (solid line), Y p = 0.41 (dashed line), Y p = 0.39 (dotted line), Y p = 0.37 
(dot-dashed line), Y p = 0.35 (double dot-dashed line) and Y p = 0.33 (three dot-dashed 
line) . 
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Figure 5: The caloric curves are shown for g2 6 Srriio4 (Y p = 0.37 - thick dashed line) and 
g2°Sni88 (Y p = 0.41 - thick full line); the excitation energies at T = 5 MeV, for different 
proton fractions (big triangles, from the left to the right Y p = 0.33, 0.35, 0.37, 0.39, 0.41 
and 0.44), the experimental data from |7]] (diamonds) and (stars), and the Fermi-gas 
law (k = 13.0 - thin dashed line) are also displayed. 
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